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■ Abstract 



A definition of surface gravity at the apparent horizon of dynamical spher- 
ically symmetric spacetimes is proposed. It is based on a unique foliation by 
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£T) • ingoing null hypersurfaces. The function parametrizing the hypersurfaces can 

o . 

be interpreted as the phase of a light wave uniformly emitted by some far- 
away static observer. The definition gives back the accepted value of surface 
5h . gravity in the static case by virtue of its nonlocal character. Although the 

definition is motivated by the behavior of outgoing null rays, it turns out that 
there is a simple connection between the generalized surface gravity, the ac- 
celeration of any radially moving observer, and the observed frequency change 
of the infalling light signal. In particular, this gives a practical and simple 
method of how any geodesic observer can determine surface gravity by mea- 
suring only the redshift of the infalling light wave. The surface gravity can 
be expressed as an integral of matter field quantities along an ingoing null 
line, which shows that it is a continuous function along the apparent horizon. 
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A formula for the area change of the apparent horizon is presented, and the 

possibility of thermodynamical interpretation is discussed. Finally, concrete 

expressions of surface gravity are given for a number of four-dimensional and 

two-dimensional dynamical black hole solutions. 
PACS numbers: 04.70.Bw,Dy 
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I. INTRODUCTION 



There has been a renewed interest in spherically symmetric spacetimes in the past half 
decade. The unexpected complexity of the problem is well illustrated by the large number 
of new 'dirty' black hole solutions. An important result proved by Visser JI] is that these 
stationary matter fields always decrease the surface gravity compared to the same mass vac- 
uum black hole. Dynamical spherically symmetric spacetimes were initially studied mainly 
to check the validity of the cosmic censorship hypothesis. More recently, supported by pow- 
erful numerical methods, considerable efforts have been focused on near-critical collapsing 
solutions at the black hole formation threshold. 

The fundamental question is what physical quantities are describing these spherically 
symmetric collapses. Undoubtedly, local quantities such as a now well defined gravita- 
tional mass function and densities belonging to the matter fields are essential. Because of 
their importance in the stationary case, we expect that generalizations of thermodynamical 
quantities will also play a major role. For stationary black holes, the surface gravity is pro- 
portional to the temperature of the Hawking radiation. On the other hand, considering the 
collapse of a spherical shell, Hiscock || proposed to identify one-quarter of the area of the 
apparent horizon as the gravitational entropy. Furthermore, Hajicek suggested that the 
Hawking effect is associated with the apparent horizon rather than the event horizon, since 
the apparent horizon in spherically symmetric spacetimes acts as the boundary of negative 
energy states. Hence we expect that some naturally generalized surface gravity for apparent 
horizons will have a crucial role as a physical quantity in dynamical spacetimes. 

It is possible to formulate a local dynamical analogue of black hole thermodynamics even 
for not spherically symmetric apparent horizons. Using a null tetrad formalism, Collins 0] 
has derived a formula for the area change of the apparent horizon, which can be interpreted 
as a generalized first law of thermodynamics. However, the temperature term in this equa- 
tion is tetrad dependent, and no unique tetrad choice is made in the nonstationary case. 
This ambiguity reflects the difficulty of selecting an appropriate distance function along the 
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apparent horizon. Hayward || uses the natural distance defined by the spacetime metric, 
although this choice is divergent in the stationary limit when the horizon becomes null. 
Hayward also presents the analogues of the zeroth and second laws of thermodynamics, and 
defines a dynamical counterpart of surface gravity, called trapping gravity. Unfortunately, 
when specializing to static spherically symmetric spacetimes again, the trapping gravity does 
not agree with the accepted value of surface gravity even for charged Reissner-Nordstrom 
black holes (see Appendix). 

There are two basic ways to introduce surface gravity in stationary spacetimes. The first, 
physically more direct method is in terms of the acceleration of stationary observers near 
the black hole horizon. This form of the definition proves to be very difficult to generalize. 
In dynamical spherically symmetric spacetimes, the observers moving on constant radius 
orbits are the most natural analogues of the static observers. However, their acceleration 
have a qualitatively different behavior, being proportional to the matter density instead of 
any possible generalization of surface gravity. 

The second, mathematically more straightforward approach is to define surface gravity as 
the inaffinity of the Killing vector field along the black hole horizon. In the general dynamical 
case the apparent horizon ceases to be null, and there are no geodesies remaining tangent to 
it. However, in spherical symmetry, the outgoing radial light rays are necessarily geodesies, 
and they are locally constant radius orbits when they cross the apparent horizon. Hence it is 
natural to attempt to generalize surface gravity as the inaffinity of these outgoing null orbits 
0. The concept of inaffinity is defined only with respect to a preferred parametrization of 
the curves. The main difficulty is how to choose this particular parametrization, considering 
that one has to get back the Killing time in the stationary limit. The normalization of the 
Killing vector is defined at spacelike infinity, which shows that our definition cannot be local 
either. 

The most important idea in this paper is to parametrize the outgoing null geodesies using 
a natural spherically symmetric foliation of ingoing null hypersurfaces. We assume that the 
labelling of these hypersurfaces is defined by the proper time of a static observer at infinity. 
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This foliation can be most easily observed by any dynamical observer, simply by observing 
a radio wave emitted uniformly by a far away reference clock. We can interpret the function 
generating this foliation as a global advanced-time cooridinate. We will see that for static 
observers in a static spherically symmetric spacetime, this advanced-time parameter agrees 
with the Killing time, (apart from a radius dependent additive constant,) which ensures that 
our generalized surface gravity indeed reduces to the accepted value in the non-dynamical 
case. 

Our natural ingoing null foliation will allow us not only to give a clear physical interpre- 
tation of surface gravity, but also to prescribe the most practical way to measure it in any 
static or dynamical spherically symmetric spacetime. Any observer moving along a timelike 
orbit can precisely measure the apparent change of frequency of a standard radio or light 
signal falling in from the far away asymptotically flat region. We can find an explicit rela- 
tion between this frequency change and the acceleration of the observer. In particular, for 
any geodesic observer crossing the horizon, the proper time derivative of the redshift of the 
infalling wave is exactly equal to the surface gravity. This is particularly interesting, since it 
means that surface gravity can be determined by performing simple frequency measurements 
only. 

In principle, there are infinitely many inequivalent ways to extend the definition of surface 
gravity to nonstationary spacetimes. One has to apply physical considerations to make the 
most appropriate choice. Actually, most of these potential extensions does not seem to 
possess any invariant physical meaning at all. In contrast, as we will see, our definition of 
generalized surface gravity is supported not only by the analysis of outgoing light rays, but 
also by simple measurements performed by both geodesic and accelerating observers. 

Another physical approach, which may lead to a different (but still not local) definition 
of surface gravity, is by using a fully dynamical generalization of the Hartle-Hawking formula 
jTj. Assuming that the apparent horizon area corresponds to the entropy of a dynamical 
black hole 0, this formula may be interpreted as a generalized first law of black hole ther- 
modynamics. It has been pointed out by Collins j4|, that the temperature term appearing 
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in this equation can correspond to some possible generalization of surface gravity only in the 
near-stationary limit. Furthermore, this temperature term can change in a noncontinuous 
way along the horizon whenever there is a jump in the matter field density. This happens 
for example at the surface of a collapsing star. In contrast, as we will see, our generalized 
surface gravity is always continuous for regular matter fields. 

Our dynamical surface gravity is defined in section H, using the inaffinity of outgoing null 
rays at the apparent horizon. We note that the same definition can be applied at any point 
of the spacetime, including the event horizon, but the physical meaning would be much less 
clear there. In section HJ, a method is described for how any observer, which crosses the 
horizon in an arbitrary way, can measure surface gravity by observing light signals falling 
in from infinity. A family of observers whose accelerations are equal to the surface gravity 
multiplied by a generalized redshift factor is also presented. In section |V|, the surface 
gravity is expressed as an integral of regular matter field quantities along an ingoing null 
curve coming from past null infinity. In section [V], an equation for the area change of the 
apparent horizon is presented. The possibility of interpreting it as a dynamical first law of 
black hole thermodynamics is discussed. In section |VT], the value of the surface gravity is 
calculated for several exact solutions. These include the charged Vaidya metric, self-similar 
scalar field solutions, 1 + 1 dimensional dilaton gravity and homogenous dust ball collapse. 
In the Appendix, while examining the properties of null congruences, Hayward's definition 
of trapping gravity is reviewed in the spherically symmetric case, and the relation to our 
formulation is discussed. We use units in which the gravitational constant and the speed of 
light satisfy G = c = 1. 

II. GENERALIZED SURFACE GRAVITY 

The surface gravity k of stationary spacetimes is defined by the behavior of the timelike 
Killing vector £ a at the event horizon. The definition has a non-local character. If £ a is a 
Killing vector field, then b^ a is also Killing for any constant b. This changes the value of the 
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surface gravity from k to bn. Therefore one must fix the normalization of £ a . The obvious 
way to do it is to require = — 1 at spacelike infinity. To calculate surface gravity, either 
one has to know the Killing vector field globally, or has to perform integration between the 
horizon and spacelike infinity to determine the 'anomalous redshift factor' [[L]]. 

There are several equivalent expressions which can be used to define surface gravity in 
stationary spacetimes. The most appropriate for generalizing into dynamical spacetimes is 

= < a ■ (i) 

This defining equation, unlike the others, only uses the value of the Killing vector £ a strictly 
on the horizon. It describes the failure of £ a being affine null geodesic. Since the wave 
vector of a light signal is affine null geodesic, k has the physical meaning of determining the 
frequency decrease, or in other words the redshift, of an outgoing light signal moving along 
the horizon, measuring it with respect to the Killing-time. Hence k describes the 'energy 
loss' of a photon trying to climb out of the black hole, but only able to move exactly along the 
constant radius horizon. No such frequency change occurs for a light signal moving exactly 
along the horizon of an extreme Reissner-Nordstrom black hole, although the redshift of a 
photon escaping from very close to the horizon to infinity can be still arbitrarily large. 

Given any spherically symmetric spacetime which is asymptotically flat at past null 
infinity, let us consider a foliation by ingoing null hypersurfaces parametrized by a function 
v. We make the parametrization of these spherically symmetric hypersurfaces unique (up 
to an additive constant) by requiring that £ a v - a = 1 at past null infinity, where £° is the 
asymptotic Killing vector. This requirement means that v is fixed by the proper time of some 
stationary observer at infinity. We can consider the function v as a global advanced-time 
coordinate. The parametrization of the null surfaces can be more conveniently fixed using 
the natural radial function p instead of £ a . At infinity £ a £ Q = —1, p ,a p ;a — 1 and £ a p ;a = 0. 
Hence in place of £ a v - a = 1 we can equivalently require p' a v- a = 1 at past null infinity. 

It is easy to see that in a static spherically symmetric spacetime £ a v - a is constant along 
the ingoing constant v lines: 
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(ev; a ) = v^v% a , p) + -C{v\ p ) a = . (2) 

Hence £ a v - a = 1 everywhere. This means that for static observers the advanced time v agrees 
with the Killing time, apart from an observer dependent additive constant determining the 
time 'zero'. 

In a dynamical spacetime the apparent horizon is not null anymore, and there are no 
geodesies which remain tangent to it. However, outgoing radial null curves are always 
geodesic because of spherical symmetry. Furthermore, since the expansion of outgoing null 
rays vanishes at the apparent horizon, the outgoing null curves are locally constant radius 
orbits when they cross the horizon. Hence instead of the Killing vector, which is very 
problematic to generalize to dynamical spacetimes, we will use the inaffinity of an outgoing 
null vector field k a to define surface gravity. The null condition and the spherical symmetry 
only fixes the direction of k a . The most difficult problem is how to fix the normalization of 
this vector field. Since we want our definition to give back the usual value for the surface 
gravity when specializing to static spacetimes, k a should agree with the Killing vector £° on 
a static horizon. We can assure this by requiring k a v- a = 1 at every point of the spacetime. 
This determines k a uniquely in a non-local way. Because k a is geodesic everywhere, the 
relation k^k a ^ = Kk a can be used to define n at every point of the spacetime which can be 
reached by an ingoing radial lightray coming from past null infinity. However, on physical 
grounds, we are interested in the value of the surface gravity only at the apparent horizon. 
Since k a p. a = only on the apparent horizon, the physical significance of k is much less 
clear elsewhere (see Fig. [1]). 

Multiplying the formula k^k a .^ = nk a by v. a , 

k = v, a k p k a . p = k 13 (v. a k a ), p - k a k p v. )a p = -k a k p v. a p . (3) 

Since v' ,l3 v. a p = v'^v.^ a = 0, for any scalar function a the vector k a = k a + av'' a will also 
satisfy k = —k a k^v- a p. This shows that the fundamental structure is not the vector field k a , 
but the function v determining the null foliation. We only have to assume that k a points in 
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a radial direction and k a v- a = 1. Since no derivative of k a appears, it is enough to choose 
any such vector at only one point, no need to construct a vector field. The vector k a can be 
not only null but also timelike or spacelike. 

Definition: Given a foliation by ingoing null hypersurfaces, parametrized by a function 
v which satisfies v- a p'' a = 1 (or ^ a v. a = 1) at past null infinity, the surface gravity at some 
point of the spacetime is defined as 

K = -k a k%, af3 , (4) 

where k a is a vector pointing in a radial direction and satisfying k a v- a = 1. 

Given any radially directed geodesic, we can parametrize it by the advanced time v. 
Then the tangent vector k a satisfies k a v ;a = 1 and the geodesic equation k^k a p = kk a . 
Multiplying by v. a , we get k = k (see Fig. [1]). 

Consequence: For any radial geodesic with tangent vector k a satisfying k a v- a = 1, the 
surface gravity k describes the inaffinity of the geodesic as 

k^k^ = nk a . (5) 

The physically most relevant case is when k a is the unique outgoing null vector crossing 
the apparent horizon and satisfying k a v- a = 1. At the horizon of static black holes this 
null vector agrees with the Killing vector, and our definition gives the standard value of 
surface gravity. The physical meaning of the dynamical k is the same as in the stationary 
case. An outgoing light signal moves along a locally constant radius orbit when it crosses 
the apparent horizon. Since the parametrization v is not affine, k determines the frequency 
decrease, that is the redshift of the light signal at the horizon (see Fig. [I]). Physically, the 
photon loses its 'energy' because of the attractivity of the black hole. 

What happens if we try to calculate the surface gravity using a different parametrization 
of the null hypersurfaces, a function v which is not asymptotically well behaving at past 
null infinity? Then we get K = —k a k l3 v. a i3 for some k a satisfying k a v. a = 1. The physical 
parametrization always can be obtained by a relabelling of the null surfaces, v = v(v). 
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Since k a v. a — Q = v', we have to rescale the vector k a and use k a = k a /v' to ensure 
that k a v. a = 1. Expressing the covariant derivative in terms of the Christoffel symbols and 
partial derivatives in some coordinate system, 

k = -k a k\ aP = —r^k a ~kP \(v'v >a ) B - TV'^l • (6) 



Hence the physical surface gravity k is related to the unphysical k as 

K = h k -7^ v " = ~, k + (zj) » ( 7 ) 



V [V 1 ) 2 V \v 

where the primes denote derivatives with respect to v. 

Given the function v, one can choose it as one of the coordinates in a null coordinate 
system x a = (v, r, 9, <fi). The metric takes the form 

ds 2 = -Fdv 2 + 2Gdvdr + p 2 dVl 2 , (8) 

where F, G and p are functions of v and r, and G > 0. The only remaining freedom is in 
choosing the r coordinate. Using the Christoffel symbols in this coordinate system, from 
we have k = ^k^T 1 ^ , and since T 1 22 = T 1 12 = 0, 

independently of the radius function p. If we choose r as an outgoing null coordinate, then 
we obtain a double-null coordinate system with F = and k = G tV /G. Another convenient 
choice is r = p, which we will use in most of the paper. 



III. PHYSICAL IMPLICATIONS 

The familiar method of determining the surface gravity of a stationary black hole is by 
measuring the acceleration of observers moving along the Killing orbits near the horizon. 
Using the coordinate system (Q) where r = p, the most natural generalization of the Killing 
vector is £ a = (1, 0, 0, 0), because it satisfies £ a p ;a = 0, and reduces to the Killing vector in 
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the static case. Since = —F, the velocity of the observers moving along these constant 
radius orbits is u a = (l/yF, 0, 0, 0). Their acceleration is a a = vPu a n, and 



In stationary spacetimes F jV = 0, and comparing with (H) we can see that this expression 
gives k 2 . In a dynamical spherically symmetric spacetime the derivative of the radius func- 
tion p = r vanishes in outgoing null directions at the points of the apparent horizon. Since 
(2G, F, 0, 0) is such a null vector field, this means that F = there. However, from the 
Einstein's equations we have F v = —87cpGT vv , and hence F <v /F and a^ap in general di- 
verges at the apparent horizon. The only combination which is always finite is (£ Q £ a ) 3 aPap. 
It is proportional to T a pC, a ^ instead of k, and zero in the static case. From these argu- 
ments we can see, that in dynamical spacetimes the acceleration of constant radius observers 
cannot be used to define any generalization of apparent horizon surface gravity. To illus- 
trate the problem more concretely, let us consider the Vaidya spacetime, for which G — 1, 
F = 1 — 2m(v)/r and p = r. Then k = m/r 2 , and 



This diverges at the apparent horizon r = 2m, whenever m! = 7^ in the dynamical 
case. 

There is a very intimate connection between the acceleration of radially moving observers, 
the advanced time v and the generalized surface gravity k. Consider an arbitrary congruence 
of curves in the constant angle radial plane, generated by some vector field u a . We do not 
assume that the congruence is geodesic, and it can be not only timelike but also spacelike 
or null. If we define the function / by / = u a v. a , the vector field k a = u a /f satisfies the 
normalization condition k a v ;Q = 1 needed in the definition of k. We define the vector field 
which is just the acceleration when u a is a normalized velocity vector. Then 
a a _ j2^./3^,a _|_ ffcafcp Multiplying this by v- a , using that the derivative of k a v- a vanishes, 
and substituting the defining relation (|J) of k, we get 




(10) 




(11) 
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a a v. a = f 2 K + u a f, a . (12) 

If the congruence is geodesic, then a a = and k = u a (l/ f)- a . This is not very surprising, 
since we have seen in the previous section that k describes the inaffinity of any geodesic 
parametrized by v. The important thing is that / has a simple physical interpretation and 
can be very easily measured. Since / = u a v ]a , the value of / gives the ratio of the global 
advanced time change Av and the observer's proper time change At along the orbit (see 
Fig 0) : 

Considering a light signal emitted by a static observer at infinity with frequency Uoo, the 
observed frequency is uj = fuj^. Introducing the redshift factor 

^= — -1 = 7-1, (14) 
u j 

the surface gravity is 

dz 

k = u a z. a = — . (15) 

This shows that for any geodesic observer, the proper-time derivative of the observed redshift 
of a standard light or radio signal is equal to the surface gravity k. Since such frequency 
changes can be very easily and most precisely determined, this is the most practical method 
of measuring surface gravity in spherically symmetric spacetimes, even in the static case. 

The proper time is measured by a clock carried by the observer. Actually, since we have 
not used that the norm of u a is —1, r does not even have to be proper time, it is enough if 
proportional to it. But to get the physical k, the frequency Uoo of the light signal must be 
determined by the proper time of a static clock at infinity. To measure the surface gravity of 
the apparent horizon, the observer must actually cross the horizon. If the apparent horizon 
is spacelike, the observer is unable to send the result of the measurement back to infinity. 
The utmost an observer far from the black hole can know is the approximate value of K 
at the event horizon, even if the physical meaning of k is not clear there. However, we 
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expect that the generalized surface gravity will play the most important role in the case of 
dynamical evaporating black hole models, when the apparent horizon is timelike and located 
outside of the event horizon. 

Equation ([P2f ) provides the most practical way of measuring surface gravity for non 
geodesic observers as well. If the norm of u a is constant along the orbits, then u a a a = 0, 



and (|12|) determines the only nonvanishing component of a a . In our null coordinate system 
let 

<=(/,^^0,o) , (16) 

assuming / > 0. Since gapu^u 1 ^ = =Fl and g a/ 3u"u+ = 0, a general observer moving in 
a radial direction can be described by the velocity u°, while is the outward pointing 
normal vector to the orbits. Defining a" = we have 



a% = 
T z + 



1-(k- u^z, p ) ul . (17) 



This follows from the fact that has to be parallel to u±, and that by contracting with 
v- a we get back (|T^). The norm of u°[_ is 1, the acceleration \a"\ can be directly measured, 
while / can be determined by observing the frequency change of light signals falling in from 
infinity. 

All observers who measure constant redshift z have acceleration proportional to the 
generalized surface gravity k. Hence to find the natural generalization of the static observers, 
one has to look for those solutions of the equation u^f- a = which reduce to the Killing 
orbits in the static limit. Unfortunately this equation is too difficult to solve in general. The 
choice / = 1/vF in (|I6| ) gives the constant radius observers studied in the beginning of 
the section. However, F is constant along the constant radius orbits only in the static case. 
Since F = on the horizon, one possible candidate for the solution would be the constant F 
orbits. Although, in general, their acceleration is not proportional to k, it is instructive to 
study these orbits in case of the Vaidya spacetime. Since F — 1 — 2m(v)/r, from u^F-a = 
we get 
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-j T = ±F T 2r-, (18) 
/4 m 



where /_ exists only for timelike, and /+ only for spacelike constant F orbits. Then 

p ( 1 \ rm" 
% \%). s = 2rm'-mF ' (19) 

At least for the linear Vaidya spacetime, where m = miv + m for some constants m 1 and 
mo, the constant F orbits have acceleration proportional to k. In the static limit, when 
m 1 — > 0, these orbits reduce to the Killing orbits. 

Although solutions of the equation uZf ;a = always exist, for more general spacetimes 
it is generally impossible find these solutions. One has the freedom to specify the value of 
/ as initial data on some surface, for example on the apparent horizon. A natural choice 
is to fix / by requiring that u% has to be tangent to the horizon. In general, the orbits 
determined by u^f- a = will not remain tangent to the horizon. Because the orbits never 
cross into the other side of the horizon, the apparent horizon will emerge as the 'envelope' 
of these orbits. Since the equation of these orbits is too difficult to solve, we can consider a 
similar equation instead, which gives essentially the same orbits very close to the horizon. 
We solve 

W%f;a = , (20) 

where 

<=(^2k'°-°) • < 21) 

and Gh is the value of G at the point of the horizon which have the same v coordinate as 
the point where the vector w^. is considered, Gh(v) = G (v,r h (v)). If we require that w^. is 
tangent to the horizon, i. e. w^F. a \h= 0, then and not only agree on the horizon, 
but even their accelerations are the same there. Defining a" = w^w^.p , and using that 
w %f-,a — 0; w2.F. a \ h — and w a d a G h \ h = w a d a G \ h , it is straightforward to check that on 
the apparent horizon 
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~a%\ h =fKwl\ h . (22) 

The orbits of would satisfy this property everywhere, but since we are interested in the 
value of k only at the apparent horizon, it is sufficient to study the much less complicated 
to" orbits. Since Gh depends only on v, the equation w"f ]a = 0, i. e. 

dvfT 2fW h drf = ° (23) 
can always be integrated. The general solution is 

r= ^w4^ iv+c(f) ' <24) 

where c(f) is an arbitrary function of the parameter / labelling the orbits. The function 
c(f) is fixed by the assumption that every orbit becomes tangent to the horizon at some 
point. For every value of v there is an orbit which is tangent to the horizon at the point 
(v, r h {v)). The value of / is set for this orbit by the condition 

dvji ( w r \ ! 

^ = WJ h = T viG h - (25) 

The function c(f) can be determined from 

1 } 1 

Th = T 2p J c~ dv + <fh) ■ ( 26 ) 



As an example, we consider Vaidya spacetime with mass m = m + m 2 v 2 . The horizon 
is at r/j = 2m, and from 



o 

,2 



1 

2fZ 



4:m 2 v = ^ • (27) 



Substituting into (BB), 



Using (|2^), we get the equation of the orbits: 



1/1 



r-2m = ± — - 2m 2 v . (29) 

2m 2 ^4/ 2 J K J 
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The vector field can be constructed by solving this equation for / and substituting into 



(pl|). For timelike horizon the orbits are always outside, while for spacelike horizon always 
inside the horizon. Each orbit becomes parallel to the horizon at a point where 1/v = 8rri2f 2 , 
depending on the value of / labelling the orbit. In the static limit when ni2 — > 0, the vector 
field u>^ becomes tangent to the constant radius orbits. 

We have seen that there always exists a family of observers instantaneously tangent to 
the horizon, such that their acceleration at the moment of touching the horizon is fhK- In 
dynamical spherically symmetric spacetimes these orbits seem to be the most appropriate 
generalization of the Killing orbits used to define surface gravity in the stationary case, fh 
can be interpreted as a generalized redshift factor. It is finite for dynamical spacetimes, but 
always diverges in the static limit. 

IV. INTEGRAL FORMULA 

In the coordinate system ds 2 = —Fdv 2 + 2Gdvdr + r 2 (dQ 2 + sin 2 8d(f 2 ), the independent 
components of the Einstein's equations are 

8vrT rr = , (30) 
rG 

2G 

8nT vr = - M r , (31) 

8irT vv = \ (FM r + GM tV ) , (32) 

/ r 2 \ 2 Vf 

8^ [Tee - ~T a a j = —K r + - , (33) 

where 

and T Q/ g is the stress tensor of matter fields. Using the radius function p, the local mass M 
can be expressed in a coordinate system invariant form as M = p(l — p' ,a p- a )/2. 
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Denning the vectors £ Q = (1, 0, 0, 0) and £ a = (0, -1, 0, 0), we have T vv = £ a £ p T a p, T vr = 
-i a ^T al3 and T„ = l a l p T a p. The vector f can be easily constructed in any coordinate 
system, since it is future directed null and £ a p- a — — 1. However, for £ a , only its direction 
is fixed locally by ^ a p- a = 0, the normalization £ Q £ Q = —F is known only after globally 
constructing the asymptotically well behaving foliation given by v. The best one can do 
locally is to define the Kodama vector || ( a by ( a p ]a = and ( a ( a = — — 1 = J?. Then 
^ a = G( a . Since the cotangent vector p. a = (0,1,0,0) can be easily constructed in any 
coordinate system, it is useful to write the more covariant forms instead of (|3lD and fl32j): 

= %~™, (35) 

^ =-^ + ^{ K -—) ■ (36) 

Using (|30D , and assuming that G approaches 1 at past null infinity, G can be expressed 
as an integral of local quantities along an ingoing radial null line: 

oo 

In G = -4tt f rT rr dr . (37) 



r 



This shows that G = 1 in the whole outer vacuum region. If the matter fields satisfy the 
weak energy condition, then G is non-increasing in ingoing null directions, and < G < 1. 
After calculating G, even F can be determined locally from the expression fl34j) of the local 
mass, F = G 2 (l-2M/r). 

From equation (|33"D, k can be expressed as an integral along an ingoing radial null curve 

oo 

« = /G(^-^r w + ^)dr. (38) 

r 

Since G is not a local quantity, before calculating k one has to evaluate the integral (p7|) to 
get G at every point of the null line. Using the expression for T vr from 

2 F 2 

T a a = "qTvt + ~Q{F rr + ^-Tee , (39) 

we can get 
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G-j) = 2G (nT rr - nT a a + ^T ee - ^) , (10) 



r J r \ T r J 

which gives another integral formula for k: 

oo 

K = G—+j2nG (T rr + T a a - -T w ) dr . (41) 

r 

This latter expression is especially useful if there are vacuum regions. 
Using equation (p5|). 



^ G (4„T/ + 5). (42) 

The best we can do for the G tV /G term in the expression (Q) of k is to take the derivative 
of (0). We obtain 

K = G^rT/ + ^)-4nJr^dr. (43) 

r 

After G is already known, the partial derivative can be written in a coordinate system 
invariant form, as a Lie- derivative along the vector field £ Q = G( a : 

dT r „ 



rr _ l)OL 



i p C G <T afi . (44) 



dv 

In the static case £° is the Killing vector, and the integral term vanishes. 

It is instructive to introduce a basis carried by observers moving along the constant 
radius orbits. Setting n a = -j^P;a and t a = -^p£, a , we have n a n a = 1, u a u a = — 1 and 
u a n a = 0. This basis is valid only outside of the horizon. The measured energy density is 
£ = Tapf^t 13 , the radial energy flow is S — Tapp^n 13 , and the radial pressure is P = T a ^n a rfi . 
We have the ingoing null vector t a — n a = ^-£ a . Then 

G 2 F O CT" _i_ p 

T rr = T a ^ = ^-T aP (t a - n Q )(^ - n p ) = ^ . (45) 
r 1 

r 

This shows that e — 25 + P must approach zero at the horizon. Similarly we can get 
T vr = G(S - e), T vv = Fe, T r r = P - S and T rr = P (l - 2M). Substituting into (0) 

lnG=-47r/ r(£ ~^ P) dr. (46) 

r r 



In the static case S = 0, and this reduces to the formula determining the 'anomalous 
redshift' <ft = —InG given by Visser From (fEj|), using that on the horizon r = 2M and 
e — 2S + P = 0, we get that the surface gravity of static black holes is |IJ 

K = £(i_ 8 7rr 2 e) . (47) 



V. AREA LAW 

Trying to obtain a dynamical analogue of the second law of thermodynamics, we calculate 
the advanced time derivative of the apparent horizon area. The radius and the local mass 
of the apparent horizon is related by = 2M^. Since 

dr h dM h ( dr h \ 

— = 2—— = 2 M v + M rT - , 48 
dv dv \ dv J 

we have 

dr h 2M, V 



dv 1 - 2M, 



Since at the horizon F = and T OT = GT r r } from (|3lD and (j32j ) we get 



The change of the horizon area is 



Using (|43f) , we get 



<fr> ""'"dv G ft (47TT&T/ + ^ 



9^ = 4vrr^ , (52) 
dv 



where 



(49) 



M r = -47rrgT/ , M„ = — ^T OT . (50) 



S7rr ^ = 7TT^ m _ , i \ • ( 51 ) 



B^,(W + ^)^(« + 4./^|. (53) 
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We can also see from (|42|), that in the null coordinate system where p = r, we have © = 
7^=7. In the quasi- stationary limit the integral term becomes negligible, and we obtain an 
expression corresponding to the Hartle-Hawking formula |7]]. If we identify one-quarter of 
the apparent horizon the gravitational entropy 0, then we may interpret equation 

(|52|) as a generalized first law of black hole thermodynamics. 

One of the problems with the temperature term is that it can be a non- continuous 
function along the apparent horizon if there is a sudden change in the matter density. 



Whenever there is a jump in T rr , the integrand in (|53|) becomes unbounded, and B stops 
being continuous too. On the other hand, since every quantity remains regular in the 
integral form (^) of k, our generalized surface gravity is always continuous when the energy 
densities are bounded. Another difficulty is that O is not necessarily positive. Since the 
radius function p is always constant in the outgoing null direction at the horizon, -p and 
hence ^r 1 is always positive for spacelike, and negative for timelike apparent horizons. 
If the weak energy condition holds, then T vv > on the right hand side of (0), and 
O > for spacelike while © < for timelike horizons. Furthermore, under the weak 
energy condition, spacelike horizons are outer, while timelike horizons are inner according 
to the classification of Hayward || (see Appendix). Hence one would expect that horizons, 
separating an asymptotically flat region from the black hole region, are always spacelike. 
However, as we will see in the next section when studying the example of pressureless dust 
collapse, near the center the apparent horizon can become timelike (see Fig. |]). This timelike 
region is separated form the outer spacelike region by points where the horizon is ingoing 
null. In black hole evaporation models the energy condition is violated, and if T vv < then 
O > for timelike apparent horizons. 

One would expect that the change of the black hole mass appears in the first law of 
thermodynamics. Instead, the right hand side of ([52]) describes the ingoing energy flux 
across the apparent horizon. Unfortunately, there is no direct relation between this energy 
flux and the change of the local mass along the horizon. Because = 2Mh always holds on 
the horizon, 

20 



8n 2rh dv dv ' 

which is independent of the surface gravity. On the other hand, the local mass at past null 
infinity can be constant even if there is some non-massless matter falling in along timelike 
orbits. 

From (|32|) we can see that the horizon value of T vv is proportional to the derivative of 
the local mass in the constant radius outgoing null direction. Hence 

= GhM)V . (55) 

dv 

Apart from the Gh correction factor, this appears to be more similar to the expected form 
of the first law, but unfortunately the derivatives in the two sides of the equation are taken 
in different spacetime directions. The derivative of the horizon area is calculated along a 
vector tangent to the horizon: 

^t = y ° Ah ' a ' 2/Q= ( 1 '^'°'°) • (56) 

In general, the norm y a y a is not constant. When the horizon is spacelike, there is a unique 
outgoing unit- vector z a tangent to the horizon. Since F = at the horizon, using ( (49| ) . 



Although this is a natural local specification of distance along the horizon, it has the disad- 
vantage of diverging in the static limit. Using (|43|), ((50|) and fl53|) we get 



^ z a A h , a = AnrJ^-T vv . (58) 



We can also write the right hand side into the form 2JirrhM )V . The fact that the tem- 
perature term appears under a square root follows from the unnatural normalization of the 
vector z a . Substituting from eq. ( |A15| ) in the Appendix, we get the form of the first law 
given by Hay ward 0]. 
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VI. EXAMPLES 



The simplest spherically symmetric dynamical spacetime for which we can calculate the 
generalized surface gravity is the charged Vaidya metric ||, describing a massless, charged 
null fluid falling into a charged black hole. In the coordinate system (0) we have G — 1, 
F = 1 — 2m(v)/r + e(v) 2 /r 2 and p = r, and it follows from (|9|) that k = m/r 2 — e 2 /r 3 . Since 
the radius of the outer and inner apparent horizons is r± = m ± \/m 2 — e 2 , the horizon 
surface gravity is 

2r± I r\ J r\ 



K± = — [l- — \= ±— Vm 2 - e 2 , (59) 



in agreement with the Reissner-Nordstrom value. This local agreement is due to the fact 
that G = 1 and the G }V /G term is the only one in the expression @ of k from where v 
derivatives could appear. We can also see that the surface gravity is always positive for the 
outer and negative for the inner horizon. Taking the partial derivatives of (0), 

8k± =Fe 3 dn_ 1 / m \ 

de r±y/m 2 - e 2 ' dm r 2 \ v / m 2 _ e 2 y 

g«+ = 4e 2 - 3m 2 

9m r 2 + ^m 2 - e 2 (m + 2y/m 2 - e 2 ) ' 

This shows that charging this type of black holes always decreases their outer-horizon surface 
gravity. If the infalling matter is not charged and satisfies the energy conditions, then 
> 0, and the inner-horizon surface gravity is always a decreasing function of time. 
The outer-horizon surface gravity also decreases for not very strongly charged black holes 
which satisfy 4e 2 < 3m 2 . However, if we interpret the solutions with ^ < as black hole 
evaporation models, the surface gravity of these not strongly charged evaporating black holes 
necessarily increases until their mass reduces to 3m 2 = 4e 2 . 

Our next example is the Roberts solution [TIJ, describing the self-similar collapse of a 



massless scalar field. The solution can be most conveniently given in the double null form 
ds 2 = —hdudv + p 2 dQ 2 , where h = 1 and p 2 = [(1 — p 2 )v 2 — 2vu + u 2 } /4. For parameter 
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values p > 1 , this solution describes the formation of an unboundedly increasing mass black 
hole, with a spacelike apparent horizon at u — (1 — p 2 )v . Since from eq. (|9]) we have 
K — h^ v /h = everywhere, the surface gravity is zero all along the apparent horizon. This 
indicates that the Roberts solution describes an extreme black hole, analogously to the 
e = m Reissner-Nordstrom metric. 

However, not all self-similar black holes has vanishing surface gravity. There is a confor- 
mally coupled scalar counterpart of the Roberts solution JTT|. The two metrics are related 
by a conformal transformation: ds 2 = —hdudv + p 2 dQ 2 , where p 2 = hp 2 , 



i 



h=-[M^+M^k) , M=- — ■ (62) 
4 V / ' u-(l-p)v K ' 



The new apparent horizon is determined by p ;V = 0. There, using (|), 

K = l k - = J (?) v = ~7 ^ ^ = 27 ^ (1 " V> ^ • (63) 
The horizon exists for p > 1, and it is a self-similarity line given by u = c(l—p 2 )v/A, where c 
is a constant weakly depending on p: 2.535 < c < 2.6667. The positivity of c shows that the 
apparent horizon is a spacelike hypersurface. Substituting into the surface gravity formula, 

8(4 -c) 

K ~ v[8(c-2) + (p^-l)c 2 } > ■ ( } 

After the moment of black hole formation, the surface gravity gradually decreases to zero 
from an infinitely big initial value, as the mass increases unboundedly. 

Our third example is the 1+1 dimensional dilaton gravity proposed by Callan, Giddings, 



Harvey and Strominger [12 1. It is defined by the action 



5 = ^/rfa; 2 v ^|e-^[ J R + 4(V0) 2 + 4A 2 ]-^(V/ J ) 2 | , (65) 

where g a p is the two-dimensional metric, is the dilaton field, ft are matter fields, and A 
is a cosmological constant. It is convenient to use a double-null coordinate system (x + ,x _ ), 
and denote the only nonvanishing component of the metric by g + _ = — \e 2p ■ (Now p does 
not denote the radius function as in other parts of the paper.) Considering \e~^ as a radius 
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function, one can interpret the solutions of this two-dimensional theory as four- dimensional 
spherically symmetric spacetimes with metric ds 2 = —e 2p dx + dx~ + j^e~ 2 ^dVL 2 . At the 
points of the apparent horizon T^e -2 ^ is locally constant in the outgoing null direction, that 
is d + cj) = 0. It follows from the field equations that <9 + <9_(p — <p) = 0, and hence by rescaling 
the coordinates, one can always introduce a coordinate system where p — (p. The existence 
of this coordinate system actually follows from the traceless nature of the two-dimensional 
stress tensor belonging to the matter fields f. t . 

The static vacuum solution of the model is given by e~ 2p = e~ 2 ^ = y — \ 2 x + x~ , where m 
is a parameter giving the mass of the black hole. The asymptotically normed Killing vector 
of this solution is £ a = (Xx + , —As"), and the horizon is at x~ = 0. Using the relation ([[]) 
defining the surface gravity of stationary black holes, we get k = A, independently of the 
black hole mass. We will see shortly, that our generalized surface gravity always agrees with 
the cosmo logical constant, even for nonvacuum dynamical solutions. The parametrization of 
the null foliation given by x + is not well behaving at past null infinity, since ^ a V a x + = Aa; + 
is not constant there. Introducing another parametrization v defined by e = Xx + , we 
have £ a V Q f = 1. Hence this new parametrization is the one appearing in our definition of 
generalized surface gravity. Suppose that we are given any asymptotically flat nonvacuum 
dynamical solution of the field equations. There always exists a coordinate system (x + , x~), 
in which p = <fi. We assume that the solution approaches the static vacuum solution at 
large distances, which means that the parametrization defined by e Xv = Xx + has to be 
asymptotically well behaving for any solution. In the (v,x~) coordinate system is the 
same as before, but the metric components are different: 



Since on the apparent horizon d + <j) = 0, that is d v g + - = 0, from ([|) the generalized surface 
gravity is 



dx + 




(66) 



9v- = 



dv 



9+ 



k = d v g 

9v- 



v- — 



A, 



(67) 
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for any dynamical solution of the theory 



There exists a semiclassical model proposed by Russo, Susskind and Thorlacius ||l3 
which reduces to the previously discussed CGHS model at the classical level. It is defined 
by the one-loop effective action 



S=}-Jdx 2 e~ 2 * (2d+d-p - 4«9+0<9_0 + \ 2 e 2p 

1 JL / N \ 

+ 2 E d+fid-fi ~ f ^ - 2j (d +P d^p + <f)d + d_p) 



(68) 

Similarly to the previous model, the existence of a coordinate system (x + ,x~) where p = 
(j) follows from the field equations. The vacuum solution has the form e~ 2p = e~ 2< ^ = 
—X 2 x + x~, again. By the same argument as in the previous paragraph, we can show that 
any asymptotically flat dynamical black hole solution of this theory has constant surface 
gravity, k = A, on the apparent horizon. 

As our last example, we calculate the apparent horizon surface gravity of a uniform 
density dust ball collapsing from the rest. The internal region is equivalent to a part of the 
Friedmann cosmology, 

ds 2 = a 2 (-drj 2 + dip 2 + sin 2 ^dft 2 ) , (69) 

where 

^\ 2 



a = c |^cos — J (70) 

and c is some constant. Since we match to a Schwarzschild solution at the worldline of 
a dust particle at ip = ipo, the coordinate values are restricted into — n < i] < tt and 
< ^ < -00 < tt/2. The proper time is r = c(r] + sinr/)/2. Since the radius function is 
p = a sin?/', the local mass is given by 

M = P -{\ -p> a p-, a ) = \ sin 3 V». (71) 

This shows that the mass parameter of the external solution is m = ~ sin 3 ipQ. If we denote 
the maximal radius of the ball by r , then r = csin^ , and 
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Since at the apparent horizon p = 2M, for 77 > the horizon is represented by the timelike 
surface rj = ir — 2ip. This timelike horizon is a future inner trapping horizon in Hayward's 
classification |5|] (see Appendix). Introducing v = r] + ip as a null coordinate, we have 
F = G = a 2 in (^), and using (^), the surface gravity belonging to this parametrization can 
be calculated as 

k — — (a^ + 2a t v) = — cot ip . (73) 

Unfortunately, the parametrization v is not asymptotically well behaving when continued 
into the vacuum region. Hence we will have to use the transformation formula (0) to get the 
physical surface gravity k, where v = t + r + 2m In — lj is the regular null coordinate in 
the external Schwarzschild region. Looking from the vacuum region, the matching boundary 
is generated by radial null geodesies with maximal radius r . The radius and proper time 
of such geodesic can be expressed using a parameter rj as 



r = r (cos^) , r = iw-^-^ + sinr/) . (74) 



2) 1 2V2m 

Since r and r has to agree at both sides of the boundary, the parameter rj also agrees with 
the inner time coordinate 77. At the boundary v = 77 + ipo, and hence 

, _ dv dv dv dr dv dt dr r o ( cos 2) 

dv dr] drdr] dt dr dt] sin^o cos (v>o - f ) 

After taking the derivative of the reciprocal of this expression, we are ready to substitute 
into (^). Here r] is the coordinate value where the constant v null line crosses the boundary. 
Since v = (71 — 2ip) + ip at the horizon and v = 77 + ipo at the matching surface, we have 
rj = 7r — ip — ipQ. The final result for the surface gravity is 

k = [ cos (^ _ ^ ) + 2 cos(2t/; ) - 3 cos(?/> + %)} . (76) 

4r sin^ (sin 

Here, the radius at the point of the horizon where k is calculated is 
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/ v^* 

r = c sin 3 ib = \ sin 3 ib , (77) 

V 2m 

and sin 2 ^ = 2m/r . At the surface of the ball, where ip = ipo, the the surface gravity 
becomes j-, which shows that our generalized surface gravity indeed changes continuously 
along the apparent horizon. Close to the central singularity, when r and ip is small, 

K ^ -1 cot ^(1 + 2 cos ^o) <0 , (78) 
the surface gravity diverges to minus infinity. 



VII. SUMMARY AND DISCUSSION 

In this paper we have proposed a generalized definition of surface gravity on the apparent 
horizon of spherically symmetric dynamical black holes. Since in stationary spacetimes the 
surface gravity is not a local quantity, our definition cannot be local either. The necessary 
non-local structure is an asymptotically regular foliation by ingoing null hypersurfaces. The 
resulting dynamical surface gravity is proportional not only to the frequency decrease of the 
outgoing light rays, but also to the acceleration of some special family of observers. Fur- 
thermore, any observer can easily measure it by observing the apparent redshift of standard 
light signals falling in from infinity. The fact that the surface gravity can be expressed as an 
integral along an ingoing null line shows that it is always a continuous function along the 
apparent horizon, even if there are sudden changes in the matter field density. 

We have also seen that the area change of the apparent horizon, which may be essential 
in possible thermodynamical interpretations, becomes directly proportional to the surface 
gravity only in the stationary limit. On the other hand, although it is well known that 
stationary black holes emit Hawking radiation with temperature proportional to their sur- 
face gravity, it is unclear whether or not a dynamical analogue of this statement can be 
formulated. There have been attempts to define dynamical temperature only at the (ap- 
proximate) event horizon of the Vaidya spacetime [|14|| . If there was discrepancy between 
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the temperature and the surface gravity, it might be linked with the non-thermal nature of 
the Hawking radiation. 

Based on the study of the examples in the previous section, we can have a number of 
conjectures on the general dynamical behavior of surface gravity in spherically symmetric 
spacetimes. It is natural to expect that the surface gravity of evaporating black holes is 
always a positive and non decreasing function of time. This case is especially important, 
since if similarly to the quasi-static limit there was a close relation between the temperature 
of the (thermal part of the) Hawking radiation and the surface gravity, then this would be 
the strongest support in favour of our definitions. Although these kinds of calculations are 
extremely difficult to perform in four- dimensional Einstein theory, it is very encouraging that 
the generalized surface gravity of black holes in the exactly solvable two-dimensional dilaton 
gravity models (CGHS and RST) is a positive mass independent constant, in accordance 
with the Hawking temperature calculations Jl5] . 



If the matter fields falling into the black hole satisfy the energy conditions, then according 
to the classification of Hayward (see Appendix), the apparent horizon is either spacelike 
outer or timelike inner. From the examples it seems very likely that spacelike outer horizons 
always have non-negative surface gravity which decreases in the outgoing direction in most 
of the physically relevant cases. The surface gravity of timelike inner horizons is probably 
always a decreasing function of time if the singularity is in the future. Since the surface 
gravity is continuous, it is initially positive even in the inner region. If this inner region is 
large enough, then the surface gravity can become negative there. 

It is possible that the positivity of k may be proved somehow by the integral formula 
(j38|), if one uses the energy conditions and that 2M < r outside of the horizon. Although 
equation (|33| ) gives us the derivative of k in the ingoing null direction, there are no similar 
relations for the other directional derivatives. Substituting (^) and fl3l|) into (|33|), we get 



SttG 4ttF 2G(M\ 
«, r = _T w -— r fy + T ( T ^ . (79) 

At the apparent horizon F = and r = 2M. Since outside of the horizon r > 2M, the 
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third term is negative for spacelike outer and positive for timelike inner horizons (see Fig 
|3]). At the boundary of these two regions, where the horizon is ingoing null, K >r is exactly 
the horizon directional derivative, and its signature is determined only by the signature 
of the angular directional pressure ^Pee- In particular, for collapsing dust Tgg = and 
K >r = 0. This indicates that in case of dust collapse the surface gravity takes its maximal 
value exactly where the apparent horizon becomes an ingoing null hypersurface. 
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APPENDIX: 

In this appendix we review some important properties of radial null congruences in 
spherically symmetric spacetimes. We mostly follow the approach of Hayward || , with the 
exception that we do not assume that our vector fields are generated by null foliations. De- 
note the future directed null vector fields generating the congruences by fc? and fc", pointing 
in the outgoing and ingoing radial null directions respectively. Define the normalization 
function / by k+k_ a = —e~*. Because of spherical symmetry, k± are geodesies, although 
they are not necessarily affinely parametrized: 

fcjjfe^ = b + k% , kZkZ. p = b^k a _ . (Al) 

Given two intersecting ingoing and outgoing foliations of null hypersurfaces determined 
by constant values of £ + and £~, there are two obvious ways to define the null vector fields. 
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The first is to set k± a = — Then b± = and fc± are affine geodesies. The another way 
is to define the null vector fields, by 

In this case the inaffinity parameters are 

b± = -kU, a , (A3) 

which can be easily checked in the null coordinate system adapted to the foliations. Note 
that in this case = — e^, in accordance with Hay ward's notation. 

The tensor 

Kf3 = 9a/3 + e f (k +a k_fs + k_ a k +l3 ) (A4) 

acts as a projection operator into the 2-spheres. For the calculation of the expansion, shear 
and twist, we need to evaluate the tensor B^p = h2 t h 5 f3 k± a;f 3. Since B^p is symmetric and its 
trace-free part vanishes, the twist and the shear is zero. The expansion is 

8± = B± a = -klp, a = -klA ]a , (A5) 
P 71 

where the area function is A = Anp 2 . Another useful formula relating the Lie derivative of 
h a/ 3 to the expansion is Q± = h al3 Ck±h a /3/2. 

Using the Einstein's equations, one can derive two useful expressions for the directional 
derivatives of the expansions. The formula corresponding to the Raychaudhuri equation is 

k±®±-,« = ~Qi + b ± Q ± - &TtT a pk a ± kl , (A6) 

where T a/3 is the stress tensor of the matter fields. The cross-focusing equation gives the 
derivative in the another null direction: 

fc»e ±iQ = -0+0- - (fc«/ iQ + b T ) 6± - ^e~ f + ^T^k%kt . (A7) 

If either of the two expansions 0+ or 0_ vanishes on a sphere of symmetry, the sphere 
is called a marginal sphere. The closure of a hypersurface foliated by marginal spheres is 
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called a trapping horizon. A marginal sphere on the horizon with + = is called future if 
©_ < and past if 0_ > 0. It is called outer if k"Q + . a < 0, and inner if k°_Q + . a > 0. If the 
weak energy condition holds, outer horizons are spacelike or null, while inner horizons are 
timelike or null. In both cases, they are null in the k" direction if and only if T^k^k^ = 0. 

The naming 'inner' for horizons satisfying k°_Q + - a > can be misleading though. These 
inner horizons can separate a trapped region from an asymptotically flat untrapped region. 
A future horizon which is a smooth connected hypersurface can be outer in one region 
and change to be inner in another region, simply by becoming timelike through ingoing 
null directions. For example in certain cases of pressureless dust collapse, the horizon can 
be timelike-inner in a region close to the regular center, analogously to the cosmological 
horizon in a collapsing universe. Going outwards, this horizon becomes ingoing-null at a 
two-sphere, and then it is spacelike-outer. Asymptotically, in the Scwarzschild region, the 
horizon becomes null again, but then in the outgoing direction (see Fig. |3[). 

Following Hayward ||, we define the trapping gravity of an outer trapping horizon by 
the formula 



Because of the factor, k h is independent of the product k°^k_ a = —e~f, giving the same 
value as if k°_ was chosen to satisfy k"k_ a = —1. However if we rescale the null vectors by 
a function a as k® = ak^_, k® = k°:/a, keeping unchanged, 9 + becomes 6 + = aO + , and 
k h turns into 



This shows that kh is invariantly defined only on the trapping horizon, where + = 0. 



pressed in a local way by the radius function p and the stress tensor of the matter fields: 




(A8) 




(A9) 



Using the cross-focusing equation (|A7|), the trapping gravity of the horizon can be ex- 




(A10) 
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In the vacuum case we get the familiar ^- value, agreeing with surface gravity of the 
Schwarzschild solution. While k,h is defined on the trapping horizon of dynamical space- 
times, the surface gravity is defined on the event horizon of stationary solutions. For sta- 
tionary spacetimes the two kind of horizons coincide. However, in general, the value of the 
trapping gravity kh is different from the value of the surface gravity. This is the case even 
for the Reissner-Nordstrom solution, the surface gravity of which is 

z r 

h *<~h \ A. 

where m is the mass, e is the charge parameter, and r/j = m + \Jm 2 — e 2 . The value of the 
trapping gravity is 

1 



(A11) 



"h = - <? ■ (A12) 

For fixed m, k is a monotonically decreasing function of e, while Ku is not monotonic. 
K — Kh only for e = and e = m. 

Working in the null coordinate system (^) where r = p, we can choose 

K = (l, ^> 0, o) , k°L = (0, -1, 0, 0) . (A13) 



'2G' 

Then e~f = G and 0+ = Since F = on the horizon, 



K » = hu F -- (A14) 

Comparing with (|42|) and (^), we obtain the relation between our generalized surface gravity 



k, the trapping gravity k^, and the temperature term B in equation ([52]): 

7 8T 

k + 4vr / r—^dr = 2r h Gn 2 H = 8tt6 , (A15) 
J ov 

rh 

where the integral is calculated along a constant v ingoing null line. We can see that the 
surface gravity agrees with the trapping gravity only in some special cases. This happens 
for example for the Vaidya spacetime, where G = 1 and T rr = 0. 
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FIGURES 
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FIG. 1. An illustration of our dynamical surface gravity definition. The function v parametriz- 
ing the null foliation is fixed by the proper time r of a far away reference clock. The velocity vector 
of the clock agrees with the asymptotic Killing vector £ a , and £ a v- a = 1. The surface gravity 
describes the inaffinity of both k a and k a (black arrows), as long as they are tangent to a geodesic 
and k a v- a = k a v- a = 1. The white arrow represents the affine geodesic wave vector l a of some 
outgoing light signal crossing the horizon. The 'frequency' l a v- a is a decreasing function of the 
advanced-time v. 
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FIG. 2. Given an observer with velocity u a , define the function / = u a v- a . If there is an 
infalling light wave emitted by a far-away static clock with frequency u^, the observer measures 
u = f^oo frequency. For a geodesic observer, the surface gravity n is equal to the derivative of the 
redshift z with respect to the proper time r, i. e. k = u a z- a where z = 1/f — 1. 




FIG. 3. Conformal diagram describing the collapse of an inhomogeneous dust ball. The 
apparent horizon is timelike inner between the points A and B, spacelike outer between B and 
C, and null in the vacuum region. If there is a sudden change in the density, for example at the 
surface of the ball, then the direction of the horizon changes non-continuously there. Depending 
on the the initial density distribution, a null singularity may appear at the center, and near it the 
apparent horizon has to become spacelike outer again. 
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